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Table 1: List of ensembles used in this study, with uq the tadpole factor and r\ja the scale from the heavy 
quark potential. The r\ja values shown come from a smooth interpolation. 

1. Introduction 

The MILC collaboration has been carrying out simulations of 2+1 flavor lattice QCD with 
an improved staggered quark action for about 10 years. The physics program has recently been 
reviewed in An important aspect of the MILC collaboration's research program has been the 
study of the light pseudoscalar meson sector. Here we give the latest update of this program. 
Compared to the last status report in [§] lattice ensembles with smaller lattice spacings, smaller 
light quark masses and lighter-than-physical strange quark masses are analyzed. Furthermore, we 
do fits based on both SU(2) and SU(3) chiral perturbation theory (ZPT), rather than just SU(3) as 
before, and we now include NNLO chiral logarithms. 

2. The ensembles and the fitting procedures 

The MILC collaboration has generated lattice configuration ensembles at six different lattice 
spacings, ranging from a w 0. 1 8 fm down to a « 0.045 fm. In the present analysis, only the a « 0.09 
fm ("fine"), a « 0.06 fm ("superfine") and a « 0.045 fm ("ultrafine") ensembles are considered. 
With our very precise numerical data, adding in coarser lattice spacings would require inclusion of 
higher order discretization effects in the fits, which is currently not feasible. 

The ensembles considered in this study are listed in Table [j]. am' is the simulation light quark 
mass, with up and down quark masses being equal, and am' s is the simulation strange quark mass. 
Notice that several ensembles have an unphysically light am' s , about 60% of the physical strange 
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Table 2: Masses (in MeV, using r\ = 0.31 17 fm for the scale) for the lightest sea-quark pions of various 
tastes at each lattice spacing. The Goldstone pion is the taste pseudoscalar and has the lightest mass of all 
tastes, while the taste singlet has the heaviest mass. The root -mean-squared (RMS) mass is the average that 
is used in the NNLO chiral logarithms. Unless otherwise indicated, the masses given are also the lightest 
valence-quark pions on each ensemble at that lattice spacing. We drop the a w 0.15 fm and a » 0.12 fm 
ensembles from the current analysis because of the large splittings and heavy singlet pions. 

quark mass, and one ensemble has three degenerate (light) quarks. These ensembles were created 
specifically to have good control over the SU(3) ZPT fits. 

We determine the scale r\ on every ensemble from the static quark potential (see [P). The 
values listed in Table p] come from a smooth interpolation. For the analysis presented here, however, 
we use a mass independent scheme, where r\ is taken from the smooth interpolation with the quark 
masses set to their physical values. This procedure avoids spurious dependence on the quark masses 
in the ZPT fits. 

Even with the use of the improved staggered (asqtad) fermions and the fairly small lattice 
spacings considered, the taste-violation lattice artifacts are significant, and need to be accounted 
for in the analysis. We do this, as in our previous studies, by using rooted staggered ZPT forms 
(rSZPT) at NLO in our chiral fits [|3| ||]. The "rooting procedure," taking the fourth root of the 
fermion determinant when generating the lattices, is used to eliminate the unwanted tastes present 
with the use of staggered fermions. As reviewed in Ref. [jl|], recent work suggests strongly that the 
procedure does indeed produce the desired theory in the continuum limit. 

As a new feature in the present analysis, our XVI fits now include the NNLO chiral logarithms 
derived by Bijnens, Danielsson and Lahde [Q, [6|, [7|]. In contrast to the NLO chiral logs, however, 
lattice artifacts are not included in the NNLO chiral logs. Instead, we use the root mean square 
average (over tastes) pion mass for the argument of the NNLO chiral logs. This is systematic at 
this order in ZPT only if chiral symmetry violations from taste-violating lattice effects are signif- 
icantly smaller than the usual chiral violations from mass terms. That begins to be true for the 
a as 0.09 fm points, and is better satisfied for the a 0.06 and 0.045 fm ensembles. It is not true 
for ensembles with a > 0.12 fm, which is why that data is omitted from the analysis. Table ||] gives 
some representative pion masses for our ensembles. 

The SU(3) chiral fits are done in two stages. The first consists of "low-mass" fits used to de- 
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Figure 1: Low-mass SU(3) chiral fits. The left plot shows the fit giving the LO and NLO LECs. The right 
plot illustrates the convergence of SU(3) ZPT fits in the continuum, with the strange quark mass fixed at 
0.6mf ys . For this test we also included NNNLO analytic terms in the fit. 



termine the LO and NLO low energy constants (LECs), namely what we call fa and B3 (at LO) and 
the Gasser-Leutwyler parameters L,- (at NLO). Here the goal is to keep only those ensembles and 
valence points where meson masses (including kaons, which have a quark of mass m' s ) are suffi- 
ciently light that SU(3) ZPT may be expected to be rapidly convergent. In addition, taste splitting 
should be small enough that omission of taste-violations from the NNLO terms (but inclusion at 
NLO) is systematic; this, for example, is another reason to drop the a k, 0.09 fm ensemble with 
am 1 = 0.00155, am' s = 0.031. After these cuts, only the three fine and one superfine ensembles 
with m' s <0.6m^ hys are included, and the valence masses are limited by m x + m y < 0.6m% hys . The 
fits are illustrated in Fig. [j] (left). To test convergence, the full set of NNNLO analytic terms may 
also be added; as shown in Fig. [I] (right), the convergence is satisfactory. Addition of such terms 
does not improve the goodness of fit, as can be seen by comparing the confidence levels (CL) of 
the two fits in Fig. |TJ While we only show data for f n in Fig. [j], the fits include data for the pseu- 
doscalar meson masses and decay constants with all combinations of valence quark masses that 
satisfy m x + m y < Q.6m$ hys . 

In the second stage, the "high-mass" SU(3) XPT fits, all ensembles listed in Table [I] are in- 
cluded with the valence masses restricted to m x + m y < 1 .2m^ hys . The LO and NLO LECs are fixed 
at the values from the low-mass fits. NNNLO and NNNNLO analytic terms are included, but not 
the corresponding logs. These terms are needed to obtain good confidence levels, and they allow us 
to interpolate around the (physical) strange quark mass. That fact that they are required indicates 
that SU(3) ZPT is not converging rapidly at these mass values, unlike the situation in the low-mass 
case. Since the LO and NLO LECs dominate the chiral extrapolation to the physical point, the 
results for decay constants and masses are insensitive to the form of these NNNLO and NNNNLO 
interpolating terms, as long as the fits are good. The high-mass fits are used to give the central 
values of the physical decay constants and other quantities involving the strange quark mass, such 
as fi, B2 and chiral condensate {uu)2, which are defined in the two-flavor chiral limit (m — > 0, m s 
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Figure 2: High-mass SU(3) chiral fits; in the left plot selected partially quenched data points are shown, 
while in the right plot only full QCD points, i.e., points with sea and valence quark masses set equal, are 
shown. 




fixed at m p s ys ). The high-mass fits are illustrated in Fig. |[ 

Finally, we also perform SU(2) ZPT fits. Here we consider the five superfine ensembles with 
m' s f» mf and the ultrafine ensemble. The fits are systematic, with rSZPT forms at NLO and 
continuum NNLO chiral logs with RMS pseudoscalar masses as arguments. The valence masses 
are restricted to m x + m y < 0.5mf hys . The fits, and the convergence, are illustrated in Fig. |3j, where 
here we show both decay constant and meson mass. 



5 



MILC results for light pseudoscalars 



Urs M. Heller 



3. Preliminary results 



In a first analysis we use, as before, a lattice scale determined from T-splittings [Jsj] which leads 
to rf ys = 0.318(7) fm [g]. With this, we obtain 



f„ = 128.0 ± 0.3 ± 2.9 MeV 
f K = 153.8 ±0.3 ±3.9 MeV 
fx/ fx = 1-201 (2) (9) . 



(3.1) 



Here, and in the following results, the first error is statistical and the second is systematic. 

Our result for f n agrees nicely with the latest PDG 2008 value, f n = 130.4 ± 0.2MeV [|(]]. 
Since f % is our most accurately determined dimensionful quantity, we can use it to determine the 
scale. This gives rf = 0.3117(6)(^3j) fm. Redoing our analysis with this more accurate scale, 
we obtain 



fK 

h 
h 



156.2 ±0.3 ±1.1 MeV 
122.8 ±0.3 ±0.5 MeV 
110.8 ±2.0 ±4.1 MeV 
A// 2 = 1.062(l)(3) 
(uu) 2 = -(279(1)(2)(4) MeV) 3 
2L 6 -U = 0.16(12) (2) 
U = 0.31(B)(4) 
L 6 = 0.23(10)(3) 
m s = 89.0(0.2)(1.6)(4.5)(0.1) MeV 
m u = 1.96(0) (6) (10) (12) MeV 
m s /m = 27.41(5)(22)(0)(4) 



f K /f„ = 1.198(2) (tl), 

B 2 = 2.87(1) (4) (14) GeV, 

B 3 = 2.39(8) (10) (12) GeV, 

Ujh = 1.172(3) (43) , 

(mm) 3 = -(245(5)(4)(4) MeV) 3 

2L 8 -L 5 = -0.48(8)(21) , 

L 5 = 1.65(12) (36) , 

L 8 = 0.58(5)(7) , 

m = 3.25(1) (7) (16) (0) MeV, 

m rf = 4.53(l)(8)(23)(12) MeV, 

m u jm d = 0.432(1) (9) (0) (39) . 



(3.2) 



Here the NLO LECs L, are in units of 10 3 , evaluated at chiral scale m r] , and the LO LECs Bj, 
quark masses and chiral condensates are in the MS scheme at 2 GeV. For the conversion from 



the bare quantities we use the two-loop renormalization factor of []1 1[]. The resulting perturbative 
error is listed as the third error in these quantities. The subscripts "2" and "3" refer to the two-flavor 
(with m s at its physical value) and three-flavor chiral limits, respectively. The quark condensates are 
related to the LO LECs by (uu) j = —fjBj/2. Quark masses, finally, have a fourth error, accounting 
for our limited knowledge of electromagnetic effects on pion and kaon masses (see [[12|] for how 
we address this). 

We note that most of our new results agree, well within errors, with our previous analysis [Q] 
and have smaller errors. Not surprisingly, an exception are the NLO LECs, some of which changed 
considerably with the inclusion of NNLO chiral logs. Similar changes have been observed in 



continuum extractions of these NLO LECs; see for example Ref. Q13], 



For the SU(3) NLO ZPT correction to the physical pion mass [14] we find 



=0.04(5) (2) 



(3.3) 
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Using one-loop conversion formulae [15] we obtain from the SU(3) NLO LECs in Eq. (3.2) the 
scale invariant SU(2) NLO LECs [[l|] 

h =3.32(64) (45) , l A = 4.03(16)(17) . (3.4) 

From our SU(2) chiral fits we obtain the preliminary result, using the scale from T-splittings, 

f„ = 128.7 ± 0.9+^ MeV , (3.5) 



to be compared to Eq. (p.l|), and using the more accurate scale from f K , 

h = 3.0(6) (±|) , h = 3.9(2) (3) , 

h = 123.7 ± 0.8+} 3 4 MeV, B 2 = 2.89(2) (+J) (14) GeV , (3.6) 

m = 3.21(3)(+j)(16)(0) MeV , (uu) 2 = -(280(2)(+j)(4) MeV) 3 , 



to be compared to Eqs. ( p.2| ) and < J3-4| ). We observe nice agreement between our SU(2) and SU(3) 
chiral fit results. This indicates that a careful use of SU(3) ZPT may be justified in the extrapolation 
of other quantities calculated on the lattice, such as Bk- 
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